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Different scenarios of the fluctuation-induced disordering of the striped phase which is formed 
at low temperatures in the triangular-lattice Ising model with the antiferromagnetic interaction of 
nearest and next-to-nearest neighbors are analyzed and compared. The dominant mechanism of the 
disordering is related to the formation of a network of domain walls, which is characterized by an 
extensive number of zero modes and has to appear via the first-order phase transition. In principle, 
this first-order transition can be preceded by a continuous one, related to the spontaneous formation 
of double domain walls and a partial restoration of the broken symmetry, but the realization of such 
a scenario requires the fulfillment of rather special relations between the coupling constants. 
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I. INTRODUCTION 

An Ising model can be denned by the Hamiltonian 

H = J2j»^i, (1) 

(iJ) 

where the fluctuating variables (spins), oj = ±1, are de- 
fined on the sites j of some regular lattice. In the standard 
version of the model 0,0,01 the coupling constants, Jy, 
are assumed to be non-zero only when i and j are the 
nearest neighbors of each other, but in a more general 
case one can suppose that they depend on the distance 
between i and j. The models belonging to this class play 
an extremely important role in the condensed matter 
physics, because they can be used for the description of a 
huge variety of systems with a two-fold degeneracy of an 
order parameter. The best known examples of such sys- 
tems are ferromagnets and antiferromagnets with strong 
easy axis anisotropy and absorbed monolayers. 

The exact solution of the Ising model on a triangular 
lattice with the interaction of only nearest neighbors was 
found in 1950 Q - In the case of the isotropic anti- 
ferromagnetic interaction it demonstrates rather unusual 
properties. Namely, the system remains disordered at ar- 
bitrarily low temperature 0, Q and at zero temperature 
is characterized by an algebraic decay of the correlation 
functions 0, Q and a finite residual entropy per site . 
The ground states of this model can be mapped @ onto 
the states of a solid-on-solid (SOS) model describing the 
fluctuations of the (111) facet of a crystal with a simple 
cubic lattice and are infinitely degenerate. 

This degeneracy is not related to symmetry and there- 
fore in a physical situation its removal by the interactions 
of more distant neighbors should be taken into account. 
If the interaction of second neighbors (characterized by 
the coupling constant J2) is included into consideration, 
for both signs of J2 the degeneracy of the ground states 
is reduced to a sixfold one pj. In terms of the SOS rep- 
resentation |6(, the ferromagnetic interaction of second 
neighbors (J2 < 0) corresponds to the positive energy 



of a step and, therefore, leads to the stabilization of the 
flat phase at low enough temperatures @- With the in- 
crease of temperature a roughening transition takes 
place, which at a higher temperature is followed by an- 
other phase transition related to the dissociation of pairs 
of dislocations 0, ^j. Both transitions belong to the 
Berezinskii-Kosterlitz-Thouless universality class. 

This scenario follows also from approximate mappings 
[Til ll^l of the considered model onto a six-state clock 
model and has been confirmed by numerous Monte-Carlo 
simulations 0, [l2L Il3| . As is typical for low-dimensional 
systems, a mean field analysis p~4| leads to the wrong 
conclusions about the character of the intermediate phase 
or the nature of phase transitions. 

The present work is devoted to the antiferromagnetic 
Ising model on triangular lattice in which the interaction 
of second neighbors is also antiferromagnetic. For brevity 
we shall call such a system a triangular-lattice Ising an- 
tiferromagnet. The structure of the ground state of the 
triangular-lattice Ising antiferromagnet with the interac- 
tion of first and second neighbors p| is shown in Fig. ^ 
Here and below we use filled and empty circles to denote 
the spins of opposite signs. 

Although this version of the Ising model have been 
also investigated by different methods pj| 0, 0, 0] , its 
properties are not as clearly understood as those of the 
model with the ferromagnetic interaction of second neigh- 
bors. In particular, Domany et al. jl5J have demonstrated 
that the formal construction of the Ginzburg-Landau 
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FIG. 1: The structure of the ground state for Ji,2 > 0. 
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functional describing the formation of the state shown 
in Fig. ^ reproduces the Hamiltonian of the Heisenberg 
model with the cubic anisotropy, which under the renor- 
malization is transformed [19) into the discrete cubic 
model [23| . However, this hardly allows one to make any 
conclusions about the properties of a triangular-lattice 
Ising antiferromagnet (i) because the cubic model al- 
lows for three different scenarios of disordering l2l l (a l- 
though sometimes only one of them is mentioned |22fl). 
and (ii) because this approach does not take into account 
the strong chiral asymmetry |23j | existing in the prob- 
lem (see Sec.[Hj|. Additionally, such a description cannot 
reproduce the well-known properties of the model with 
only first-neighbor interaction. The same can be also said 
about the mean- field approach of Kaburagi et al. 
which gives an unrealistic prediction of a three-sublattice 
intermediate phase. 

A more direct analysis of the fluctuations which can be 
responsible for the disordering of the striped phase has 
been undertaken by Hemmer et al. |l7l llSj. However, 
these authors have considered the formation of only one 
type of domain walls, whose appearance (as is discussed 
in Sec. IHIfl can lead only to a partial restoration of the 
broken symmetry. Thus, all conclusions of Ref. ^3 an d 
Ref. are based on the essentially incomplete physical 
picture and have to be reconsidered. 

This article is devoted to comparison of different mech- 
anisms and different scenarios of disordering which are 
possible in the striped phase of a triangular-lattice Ising 
antiferromagnet. Although a number of numerical works 
[24l I2H l2t| give evidence on the existence of a single 
first-order transition (for this or that choice of the rela- 
tion between the coupling constants), in a situation when 
different mechanisms of disordering compete with each 
other one hardly can be sure about the universality of 
this result. 

Our analysis assumes that J\ is much larger then all 
other coupling constants, which allows us to use the an- 
alytical methods based on the separation of different en- 
ergy scales. This limit is also of interest because it corre- 
sponds to the case of strong chiral asymmetry and is the 
most relevant one for many physical realizations of the 
model. The outlook of the article is as follows. 

In Sec.[n]we discuss the symmetry of the ground states 
and the structure of domain walls and estimate the tem- 
perature at which the free energy of a single domain wall 
vanishes as a result of thermal fluctuations of this wall. 
In Sec. IIIII we show that the free energy of a double do- 
main wall can be expected to vanish at much lower tem- 
perature than that of a single wall, and argue that this 
suggests a possibility of a two-transition scenario. 

In Sec. IIVI the spontaneous formation of a network of 
single domain walls is analyzed. In the limit when this 
network has to be diluted, it becomes clear that it has 
to appear via a first-order phase transition. When only 
the interaction of up to third neighbors is taken into ac- 
count, the estimates for the temperatures of this transi- 
tion and of the spontaneous formation of double domain 



walls coincide with each other, which suggests that there 
is only one phase transition in the system. In Sec. [3 a 
phenomenological free-energy functional is constructed 
which allows one to confirm that the transition to the 
disordered phase has always to be of the first order. 

The interplay between the two main mechanisms of the 
disordering is analyzed in Sec. I VII whereas the results are 
summarized in Sec. I Vlll The short Appendix is devoted 
to a formal derivation of an exact upper boundary for 
the temperature at which the long-range order in Oj is 
destroyed by thermal fluctuations when J\ ~ 00. 



II. GROUND STATES AND DOMAIN WALLS 

The structure of the ground state of the triangular- 
lattice Ising model with the antiferromagnetic interac- 
tion of nearest and second neighbors is shown in Fig. 
^ whereas Fig. [^illustrates the classification of neighbors 
on a triangular lattice according to their distance from a 
given cite (denoted by zero). In the following we assume 
that J3, the coupling constant describing the interaction 
of third neighbors, can also be non-zero, but satisfies the 
constraint 27, 28], J3 < J2/2, which is required for the 
stability of the striped ground state of Fig. ^ The role 
of the interactions of more distant neighbors will be dis- 
cussed in Sec. IVII To make a situation more transparent 
we assume that the interaction of nearest neighbors, J\, 
is much larger then all other coupling constants. 

The sixfold degeneracy of the state shown in Fig. 2] 
corresponds to the violation of the Z2 x Z3 symmetry, 
where Z-j is related to the possibility of interchanging 
positive and negative spins, and Z% to three possible ori- 
entations of the stripes formed by spins of the same sign. 
Each of the six ground states can be associated with a 
unit vector pointing either in positive or negative direc- 
tion along one of the three axes, which are perpendicular 
to each other. These six directions can also be put into 
correspondence with the six faces of a cube. In such a 
representation the group Z2 is related to the reflection 
symmetry which transforms the opposite faces of a cube 
into each other, whereas the group Z3 corresponds to the 
cyclic permutations of the three axes. 

In systems with a discrete degeneracy the destruction 
of a long-range order has to be driven by thermal acti- 
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FIG. 2: The classification of neighbors on triangular lattice. 
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vation of infinite domain walls. The appearance of a se- 
quence of more or less parallel walls is expected when the 
intersections of walls with different orientations are ener- 
getically unfavored [2!j. The alternative option consists 
on the appearance of a network of intersecting domain 
walls with different orientations, which is would be fa- 
vored by a negative energy of domain wall intersections 

Fig-EI a ) shows an example of a lowest-energy domain 
wall separating two different ground states. Each segment 
of such a wall connects two neighboring spins of the op- 
posite signs and separates two second neighbors of the 
same sign. The energy of this wall per unit length, 



2 J 2 - 4 J 3 > , 



does not depend on J\ , because its presence does not lead 
to the violation of the constraint 



(2) 



on any triangular plaquette. This condition is satisfied 
when a plaquette contains spins of both signs. 

When crossing a domain wall the direction of the 
stripes formed by the spins of the same sign changes 
by 60° and, therefore, the direction of a lowest-energy 
wall is uniquely determined by the pair of ground states 
which it separates. In other terms, the energy of a do- 
main wall in the considered model is strongly dependent 
both on its orientation and on which states it separates. 
From the analysis of systems with a threefold degeneracy 
it is known that such a property (the chiral asymm etry 
jl^l) can lead to the chan ge o f the universality class j^Ij 
or even of the order [^J of a phase transition. 

This is the reason why the analogy with the cubic model 
[20II2H without chiral assymetry (which follows from the 
Landau-Ginzburg analysis of Domany et al. is insuf- 
ficient for understanding the properties of a triangular- 
lattice Ising antiferromagnet with J2.3 <C J\. 
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FIG. 3: Low energy domain walls: (a) a straight one; (b) with 
a kink. 



Any fluctuations of domain wall are impossible with- 
out the violation of the constraint J2J, and therefore re- 
quire the energies proportional to J\ ^> E dw . FigureOHb) 
shows an example of the simplest elementary defect (a 
kink) which can be formed on a straight domain wall. 
The energy of such a defect, 

E k = 2Ji - 4J 3 , 

has to include a contribution proportional to J± , because 
the formation of a kink requires to have one plaquette at 
which all three spins are of the same sign, and therefore 



= 3. 



In Fig. Gib) this plaquette is shown by bold lines. 

At finite temperatures the free energy of a domain wall 
(per unit length) can be estimated as the difference be- 
tween its energy and the entropic term related to the 
formation of kinks 30] . For T <C E k this gives 



F dw (T) » E dw - 2Texp(-£ k /T) 



(3) 



It is well known that when the intersections of domain 
walls are unfavored, one can expect the formation of a 
dilute sequence of parallel (on the average) walls when 
the free energy of a single wall, F dw (T), becomes equal 
to zero ^Em. For F dw (T) defined by Eq. © this takes 
place at 



T = Ti 



ln{E k /E dv 



(4) 



where we have taken into account that E dw <C E k (as 
a consequence of J2 <C J\). Eq. (@J shows that for 
^2,^/3 <C Ji the temperature T\ only weakly depends on 
J2 and J3 and is determined mainly by J\. 

In the limit J\ — > 00 the fluctuation-induced vanish- 
ing of the free energy of a single domain wall becomes 
impossible. This manifests itself in the divergence of the 
expression for T\ given by Eq. (0J for E k — > 00. Quite re- 
markably, even in this limit there still remain possibilities 
for a fluctuation-induced destruction of the long-range 
order. They are related to the spontaneous formation of 
double domain walls (which is discussed in Sec. IIII|I and 
of a domain- wall network (discussed in Sec. HVfl . 



III. SPONTANEOUS FORMATION OF DOUBLE 
DOMAIN WALLS 

A double domain wall consists of two parallel single 
walls [see Fig. 0{a)] and separates two ground states 
with the same direction of spin stripes. The energy of 
a straight double wall per segment, E ddw , is given simply 
by the energy of two single domain walls from which it 
consists, -Eddw = 2E dw . The interaction between two par- 
allel single walls does not appear even if one takes into 
account the interactions of spins with their fourth and 
fifth neighbors (see Fig. However, it turns out that 
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the fluctuations of a double wall cost less energy than 
those of a single wall, as a consequence of which its free 
energy can easily become smaller than that of a single 
wall. 

Fig. E[b) demonstrates that the fluctuations of a dou- 
ble domain wall are possible without the violation of the 
constraint @. From this figure it is clear that each cor- 
ner on a double wall requires the appearance of an ad- 
ditional segment of a single domain wall. The energy of 
such a defect, E c = 2J 2 , does not depend on J\ and J 3 . 
At finite temperatures the expression for the free energy 
of a double wall which takes into account t he p resence of 
thermally activated corners is of the form |18l l36| 



F ddw (T) = 2E dw - Tln[l + exp(-£ c /T)] 



(5) 



The spontaneous appearance of a diluted sequence of 
such walls can be expected to take place when F ddw (T) 
becomes equal to zero. The condition F ddw (T) = can 
be rewritten as 



T = 



■ ln[exp(2i; dw /T) - 1] 



(6) 



Note that Eq. © and, therefore its solution, T 2 , do not 
depend on J\. 

In the case of J3 = the solution of Eq. © gives 



T2 = I2J2 



where \m Li 



72 




(7) 



7.112 



whereas for E dw <C E c an expansion of the exponent in 
the right-hand side of Eq. (JSJ allows one to find that 



T 2 



E r 



2J 2 



In (E c /2E dw ) In ( J 2 /2£ dw ) 



(8) 
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FIG. 4: Double domain walls: (a) straight; (b) with two cor- 
ners. 



Comparison of Eq. (0 and Eq. (JSJ with Eq. (@J demon- 
strates that for J 2 <C J\ the destruction of the long-range 
order in oj cannot be driven by the spontaneous forma- 
tion of a sequence of single domain walls, because the 
analogous sequence of double domain walls can be ex- 
pected to appear at much lower temperature. A numeri- 
cal calculation of T 2 for an arbitrary relation between J\ 
and J 2 (at J3 = 0) in terms of the one-dimensional SOS 
model which takes into account also more complex fluc- 
tuations of a double domain wall can be found in Ref. ITU 

In the limit J\ — > 00 at low temperatures the system 
has to be completely frozen in one of its ground states, 
because any finite size fluctuation on the background of 
the striped ground state requires the violation of the con- 
straint ||2J on at least two plaquettes (see Fig.[5J). Accord- 
ingly, the exact expression for the free energy of a double 
domain wall cannot contain any additional contributions 
related to the suppression of finite size fluctuations. This 
has allowed Shi and Wortis |36] to conjecture that in the 
limit J 1 — > 00 the condition F ddw (T) — [with Fddw 
given by Eq. JSJ] determines the exact value of the tran- 
sition temperature. However, this conclusion can be valid 
only if the spontaneous formation of a diluted sequence 
of double domain walls is not preceded by the sponta- 
neous formation of a network of single domain walls (see 
Sec.lFvl. 

It is rather evident that the average direction of a fluc- 
tuating double wall containing thermally activated cor- 
ners will be perpendicular to the direction of spin stripes 
on both its sides. The spontaneous formation of a se- 
quence of such walls leads to the restoration of the Z 2 
symmetry between the two states with the same direc- 
tion of stripes and the reduction of the broken symmetry 
to Z 3 . 

Since the concentration of walls, v{T), is restricted by 
their collisions (which are responsible for the reduction 
of the entropy of their fluctuations) , Shi and Wortis [3(| 
have concluded that this phase transition is continuous 
and belongs to the Pokrovsky-Talapov |3l| universality 
class. Accordingly, in the vicinity of T 2 one should have 
v{T) cx (T - T 2 ) 1/2 . Note that the value of the correla- 
tion length describing the decay of the correlation func- 
tion (<7i0j) in the direction along the stripes is inversely 
proportional to v(T). 

On the other hand, Einevoll and Hemmer |18| have 
argued that this phase transition cannot be continuous, 
because the temperature at which F ddvr (T) vanishes is 
different for different directions of a double wall. In our 
opinion this conclusion is completely unjustified. The de- 
pendence of Fddw on the direction of the wall manifests 
itself in the spontaneous formation of a sequence of walls 
with the same (on the average) direction, and is a nec- 
essary condition for the applicability of the Pokrovsky- 
Talapov theory [3l| rather than an obstacle for its valid- 
ity. 

If one assumes that the value of J\ is finite, but large 
in comparison with J 2 , the fluctuations in the low tem- 
perature phase become possible. Each point where a pair 
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FIG. 5: Closed loops formed by (a) single domain walls, (b) 
double domain walls. 



of single walls is created (or annihilates) cost the energy 
close to 2Ji, so at low temperatures there will exist a fi- 
nite concentration of highly anisotropic loops formed by 
such walls [see Fig. Efa)]. However, at T w T 2 -C Ji the 
average distance between them will be much smaller then 
their size, as a consequence of which their presence can be 
neglected. On the other hand, the size of the closed loops 
formed by double domain walls [see Fig. [Sfb)] diverges 
when T — > T%- From the theory of the commensurate- 
incommensurate transition it is known |37| that this is 
accompanied by the change of the type of the phase tran- 
sition from the Pokrovsky-Talapov universality class [3lJ 
to that of the Ising model. However, the behavior will be 
changed only in a narrow region around T2, which will 
be exponentially small in J1/T2. 

Since the considered phase transition is related to the 
restoration of Zi symmetry, the emergence of the Ising 
critical behavior looks rather natural. The different uni- 
versality class in the case of J\ = 00 can be explained by 
the extremely anisotropic nature of domain walls in that 
limit, which prevents the merging of different domains of 
the same state. 

If the spontaneous formation of a sequence of double 
domain walls indeed takes place as a separate phase tran- 
sition it has to be followed (with a further increase of 
temperature) by a second phase transition related to the 
restoration of Z% symmetry. The completely disordered 
phase above this transition will look like a mixture of fi- 
nite domains of all six ground states. Since double walls 
do not change the orientation of stripes, the second phase 
transition requires the formation of single domain walls 
of all possible orientations. Above we have found that 
the spontaneous formation of double domain walls takes 
place when the free energy of a single wall is still much 
larger then temperature, so it looks rather plausible that 
the two phase transition may be well separated from each 
other. However, to check if it is really so, the formation 
of a network of single domain walls has to be studied in 
more detail. 



IV. SPONTANEOUS FORMATION OF A 
DOMAIN- WALL NETWORK 

Like in the previous section it will be convenient to 
start by considering the case of J\ = 00. In this limit 
all single domain walls have to be straight due to the 
absence of kinks. 

Fig. shows how such walls can intersect or merge 
with each other without violating the constraint • The 
energy of these intersections does not depend on J\ or J2 . 
In particular, the energy of the 120° junction shown in 
Fig- EI a) is simply equal to zero, E a = 0, whereas for the 
60° junction shown in Fig. EJb) it is given by = 4J 3 . 
The intersection shown in Fig.[S{c) can be considered as 
an overlap of two 60° junctions. The energy of this object 
is equal to = 12J 3 . 

An important feature (already mentioned in Sec. ITTfl . 
which plays a crucial role in determining the structure 
of a domain- wall network for J\ = 00, is that the direc- 
tion of each wall is uniquely determined by the states 
which it separates. A possible structure of a network 
which is formed by straight walls and satisfies this cri- 
terion is schematically shown in Fig. Here the letters 
A, B and C are used to denote the domains with three 
different orientations of stripes. Note that all walls be- 
tween A and B are parallel to each other. The same is 
true for all walls between B and C, as well as for all walls 
between C and A. 

For the sake of clearness we have not shown in Fig. [7| 
which of the two versions of A, of B , or of C (related to the 
change of sign of all spins) is realized in each particular 
domain. This depends on the exact positions of domain 
walls. 

The structure of the network shown in Fig. 0has been 
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FIG. 6: Three low-energy domain walls with different orien- 
tations can merge or intersect with each other. 
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chosen to maximize the entropy for the given total length 
of the walls. The network of such a kind is characterized 
by a large number of zero modes which do not change its 
energy. For example, each domain of the type A can be 
moved to the left or to the right. This changes the areas of 
all domains of the types B and C which are adjacent to it, 
but the total length of domain walls (and, therefore, the 
total energy of the network) is conserved. Analogously, 
all domains of the types B and C can be moved in the 
two other directions. When a domain is moved by one 
lattice unit in such a way, the signs of all spins inside it 
are reversed. 

A combination of all three types of zero modes allows 
to change the size of a three-domain complex (a bubble) 
formed by neighboring domains of three different types 
without changing its position. An example of such a bub- 
ble is shown in Fig.0by the bold line. The zero modes of 
this particular type are called the breathing modes [35| . 

The existence of breathing modes has been discovered 
by Villain |33j| when studying the formation of a honey- 
comb network in which each domain has the shape of a 
hexagon. In such a network the size of each hexagon can 
be changed without changing the total length of domain 
walls. A honeycomb network is formed in a system with 
a threefold degeneracy in which a domain wall of a given 
type (for example, a wall between A and B) can have only 
three particular orientations out of six that arc generally 
allowed. In these terms, in our problem a wall of each 
type can have only two orientations out of six, whereas 
in the three-state Potts model on a triangular lattice all 
six orientations are allowed for domain walls of any type. 

The entropy which can be associated with the existence 
of zero modes can be estimated as InM per mode, where 
M is the typical number of the configurations which can 
be spanned by a given zero mode. It is clear that in a di- 
luted network M has to be proportional to L, the typical 
distance between the centers of neighboring bubbles (in 
lattice units). This allows one to estimate the free energy 




FIG. 7: A possible structure of a domain-wall network 



(per unit area) of a network shown in Fig. [7] as 



Fnw(L) 



where 



V3 



(9) 



E„ 



En, + Eh 



4J 3 



The first term in the right-hand side of Eq. describes 
the energy of domain walls, the second term - the energy 
of merging points and the third term is related to the 
entropy of zero modes. The expression @ has the same 
structure as the one proposed by Villain [33l l8^ | for a 
honeycomb network. 

The variation of F nw (L) with respect to L reveals that 
this function has two minima, one of which is situated at 
L = oo and corresponds to the absence of any network 
and another at L = Lq < oo. The free energies of these 
two minima become equal to each other when 



Lo 



3T 



In L = 1 



E„ 



T 



(10) 



This shows that the applicability of Eq. (O, which as- 
sumes L> 1, requires to have -Edw <C T -C E mp , that 
is i?dw -C J-z- With the increase of £W/J 2 the value of 
Lq at the transition is decreased and for E^w ~ J2 it 
becomes comparable with 1, which means that Eq. (0 is 
no longer applicable. 

In the limit of E^ w <C J 2 (when E mp w 2 J 2 ) the tem- 
perature of the first-order phase transition related to the 
formation of a domain-wall network, which follows from 
Eqs. l(TU)l . can be estimated as 



mp 



2J 2 



\n(E mp /E, 



dw J 



ln(J 2 /E dv 



(11) 



whereas the value of Lq (which determines the correlation 
radius for the fluctuations of ctj) at the transition point 
is given by 



3 J 2 /£(Iv 



2 ln( J 2 /£ dw ) 



(12) 



With the decrease of J3 the ratio -E mp /-Edw is de- 
creased, which leads to the decrease of L c . For J 3 ~ J 2 
(that is Edw ~ E m p) the value of L c following from 
Eqs. (|10f> becomes comparable with 1, which means that 
the approach based on the minimization of F nw (L) is 
no longer applicable. However, since the decrease of J3 
makes the first-order nature of the transition more and 
more pronounced, one can expect that it will remain of 
the first order even when the formation of a domain-wall 
network does not allow for a quantitative description. In 
the next section this conclusion is confirmed with the 
help of a phcnomcnological analysis which does not take 
into account any details of a domain- wall network struc- 
ture and, therefore, is applicable in a wide interval of the 
values of J3 (including J 3 = 0). 
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The finiteness of Ji cannot be expected to be of any im- 
portance for the phase transition related with the spon- 
taneous formation of a domain-wall network. It allows 
for fluctuations of single domain walls, which no longer 
have to be straight, but, as has been shown by Villain 
|34j for a hexagonal network, this does not lead to any 
qualitative changes. 

V. HIGH-TEMPERATURE PHENOMENOLOGY 

Like above, it will be convenient to start the analysis by 
considering the case of Ji = oo . In this limit the manifold 
of the allowed states coincides with the manifold of the 
ground states of the system with only nearest neighbor 
interaction, which can be put into correspondence with 
the states of the (111) facet of a crystal with a simple 
cubic lattice @ . 

In particular, the ground states whose structure is 
shown in Fig. map onto the states with the maximal 
possible slope, corresponding to the transformation of the 
(111) facet into one of the three facets of the (100) family. 
The direction of the slope is determined by the direction 
of stripes in spin representation. These states do not al- 
low for the formation of any finite size defects, because 
this would require locally the further increase of a slope. 

In the limit of T — * oo all terms in the partition func- 
tion become equal to each other, therefore in this limit 
all correlation functions in the system with J\ = oo have 
exactly the same form as in the model with the interac- 
tion of only nearest neighbors at zero temperature. The 
corresponding phase is characterized by the zero slope 
and a logarithmical divergence of fluctuations @ of the 
discrete variable n describing the position of a surface. 
One can expect that the same phase will be also stable 
at large but finite T. 

In this phase the large-scale fluctuations of n can be 
described by a continuous free energy functional, 

F cS {n} = J d 2 vf 2 {n} , / 2 {n} = f (Vn) 2 , (13) 

in which the discreteness of n is neglected . At T = oo 
the dimensionless effective rigidity K (which is of en- 
tropic origin) is equal to K = tt/9. In terms of the SOS 
representation the energy of a step (per unit length) is 
equal to — 2J 2 , therefore for J 2 < the decrease of tem- 
perature is accompanied by a monotonic growth of K § , 
which at K = tt/2 leads to the phase transition to the 
smooth phase. 

We are now considering the opposite case of J 2 > 0, 
when the decrease of T from T = oo should be accompa- 
nied by the decrease of K from K = Kq. Since we know 
that at lower temperatures the triply degenerate phase 
with a finite slope has to be formed, the phenomenologi- 
cal functional ljT3|) has to be replaced by a more complex 
one, 



where the second term in the integrand, 

h{n} = -Jjr3[(eiV)n][(eaV)n][(e 3 V)n] , (15) 

favors a finite slope in one of the three equivalent direc- 
tions set by the three unit vectors e Q (where a = 1, 2, 3) 
forming the angles of 120° with each other. The last term 
in the integrand, fi{n], is required to stabilize a finite 
value of a slope. On general grounds, one can expect that 
the expansion of fi{n\ in powers of Vn starts from the 
fourth-order contribution: 

/4{«} = y(Vn) 4 + ... . (16) 

Note that the free energy functional l|14l) has nothing 
in common with the Ginzburg-Landau functional con- 
structed by Domany et al. |l5j . In particular, in the latter 
the third-order term can appear only in the absence of the 
particle-hole symmetry (which in terms of the spin rep- 
resentation corresponds to sj => — Sj). In the model that 
we consider this symmetry is, naturally, always present, 
but F e g {n} has to contain the third-order term just as a 
consequence of the symmetry of the problem in terms of 
the SOS representation. 

From the form of F c s{n} it is clear that a phase tran- 
sition between the phases with zero and finite slopes can- 
not occur in a continuous way. The three equivalent aux- 
iliary minima of F c g{n} are formed at finite values of 
|Vn|, and a first-order phase transition takes place when 
the decrease of K makes the free energy in these minima 
lower than in the central minimum at |Vn| = 0. 

Thus we have obtained an additional confirmation of 
the conclusion that a phase transition from the disor- 
dered phase with a zero slope to a phase with a finite 
slope has to be of the first order. Naturally, the construc- 
tion of a phenomenological functional does not allow one 
to distinguish whether it has to be a direct transition 
to the completely frozen phase with the maximal possi- 
ble slope, or a transition to an intermediate phase with a 
smaller slope (that is, with a finite concentration of spon- 
taneously formed double domain walls), which at lower 
temperatures will be followed by the second phase tran- 
sition. 

At finite J\ the height variable n can no longer be 
uniquely defined. In terms of n each plaquette at which 
the condition J2J is violated corresponds to a screw dis- 
location on going around which n changes by ±6 . The 
core energy of such dislocations is close to 2Ji, whereas 
their logarithmic interaction is too weak to keep them 
bound in pairs [8|- In such a situation the effective free 
energy should be a functional not of a multivalued vari- 
able n, but of its derivatives, 

m a = (e Q V)n, 

which in the presence of free dislocations do not have to 
satisfy the condition 



F eS {n} = / d 2 r(f 2 {n} + f 3 {n} + / 4 {n}) , (14) 



mi + m 2 + iris = . 



(17) 
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This can be taken into account by making in Eqs. 
((T3)l - ljro)) a replacement, 

2 

(Vii) 2 — [m\ + m\ + m|) , (e a V)n m a , 
and adding to it a new contribution, 

JD = ~( m l + m 2 + m 3j 

(where \nK^ cx 2Ji/T), which controls the fluctuations 
of the difference between the densities of positive and 
negative dislocations. However, the minimums of this new 
functional, which instead of the two variables encoded in 
Vn depends on the three variables m a , will be achieved 
when they satisfy the condition Q17[). and therefore the 
conclusion on the first order of the transition (obtained 
at J i = oo) will hold also at large but finite J\. 

VI. COMPARISON OF DIFFERENT 
SCENARIOS OF DISORDERING 

In the two previous sections we have demonstrated that 
the phase transition to the completely disordered phase 
(which can be associated with the formation of a net- 
work of single domain walls) has to be of the first order. 
This still leaves possibilities for two different scenarios of 
disordering. The formation of a network of domain walls 
can either happen at T — T nw > T2 as a separate phase 
transition, or at T nw < Ti. In the latter case at T > T nw 
the system is already in the disordered phase (in which 
the domains of all six ground states are intermixed with 
each other) and nothing special can be expected to hap- 
pen at T — T2. In that case the only phase transition 
takes place at T = T nw . 

In the Appendix we demonstrate how at J\ — 00 one 
can construct T+, an exact upper boundary for T c , the 
temperature of a phase transition from the completely 
frozen phase, T c = min{T2, T nw }. The expression we ob- 
tain, Eq. (|A2() , is applicable for an arbitrary relation be- 
tween J 2 and J3. In the case J3 = it gives T+ « 7.54 J, 
which is only 6% above the value of T 2 ( J3 = 0) = 72T2 
discussed in Sec. IIIII Since there are no reasons for this 
boundary from above to give an extremely accurate es- 
timate of T c , one can expect that the real value of the 
temperature at which the fluctuations on the background 
of the completely frozen phase do appear will be even 
lower than T2, which means that for J3 = the single- 
transition scenario is realized. 

Comparison of Eq. (JSJ) with Eq. (|llf> shows that in the 
region of parameters in which one can rather accurately 
estimate both T2 and T nw , these two temperatures co- 
incide with each other with the logarithmical accuracy. 
This gives a hint that both mechanisms may be different 
manifestations of the same phenomenon. 

This idea looks even more plausible when one notices 
that the fluctuations of a double domain wall considered 
in Sec. IIIII can also be discussed in terms of zero modes. 



That is, the parallelogram in the middle of Fig. 0Jb) also 
can be considered as a domain which can be moved along 
the direction of domain walls that are adjacent to it with- 
out changing the total length of the walls. Apparently, a 
typical network will have a less regular structure than 
the network shown in Fig. \7\ and will incorporate some 
fragments looking like finite pieces of a double wall. How- 
ever, the contribution to the free energy from each zero 
mode has to be of the form 

cE dw L + (E a + E h )-T hi L, 

where c ~ 1, which after summation over all modes will 
reproduce the general structure of Eq. JHJl . This suggests 
that there should be only one phase transition which is 
of the first order and is related to the appearance of a 
network formed both by single and double domain walls. 

Comparison of Eq. JSJ with Eq. shows that the 
realization of the two-transition scenario with T2 < T nw 
requires E c <C E mp . When only J\, J2 and J3 are as- 
sumed to be non-zero, E c satisfies the relation 

E c — E a + Eb + Edw > E mp , 

and, therefore, any prerequisites for the phase transition 
splitting are absent. Nonetheless, they may appear when 
one takes into account the interaction of more distant 
neighbors. 

In particular, if one includes into consideration the in- 
teraction of fifth neighbors (see Fig. it turns out that 
Edw, Emw and E c remain unchanged, and therefore T 2 
should not depend on J 5 . On the other hand, E mp is 
increased by 8J5, which according to Eq. will shift 
the value of T nw upwards. Thus, one can expect that the 
increase of J5 will lead to the splitting of the phase tran- 
sition into two. 

However, the positiveness of J4 works in the opposite 
direction and therefore in a realistic system in which the 
coupling constants continuously depend on the distance 
between sites, the influence of J5 is likely to be compen- 
sated by the influence of J4 (unless J4 is, for some rea- 
sons, of the opposite sign). Nonetheless, the possibility 
of phase transition splitting is not entirely prohibited, 
although its realization in some physical system is not 
very probable. 

VII. CONCLUSION 

In the present work we have investigated the scenarios 
of disordering of the striped phase which is formed in a 
triangular-lattice Ising model with the antiferromagnetic 
interaction of nearest and next-to-nearest neighbors. Our 
analysis has shown that the destruction of such an order- 
ing has to take place via a single first-order phase tran- 
sition. 

The nature of this transition becomes more transpar- 
ent in the case J2 — 2J3 <C J2, when it can be discussed 
in terms of the formation of a diluted network of domain 
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walls, which is characterized by an extensive number of 
zero modes. In this limit one can find how the transi- 
tion temperature and the value of the correlation radius 
at the transition point depend on the parameters of the 
model [see Eq. and Eq. I|12l) ]. In the opposite limit 
of J3 — > the transition temperature has to be propor- 
tional to J2, the only energy scale which is relevant for 
J2 <C Ji , whereas the correlation radius at the transition 
point has to be comparable with 1. 

We also have shown that the formation of a domain- 
wall network could be anticipated by a transition to 
the intermediate phase, in which the long-range order 
in terms of oj is destroyed by the spontaneous forma- 
tion of a sequence of parallel (on the average) double 
domain walls, whereas the long-range order in the ori- 
entation of stripes formed by the spins of the same sign 
still exists. This transition would be characterized by a 
combination of the Ising (in the very narrow vicinity of 
the transition temperature) and the Pokrovsky-Talapov 
(in a more wide temperature interval) critical behaviors. 
However, in a system with only three coupling constants 
(Ji, J2 and J3) no prerequisites for the realization of such 
a scenario can be found. 

Nonetheless, they may appear when the interaction of 
more distant neighbors is taken into account. In partic- 
ular, the splitting of the phase transition into two is fa- 
vored by the increase of J5. However, the possibility of 
the realization of such a scenario depends on the fine 
interplay between different coupling constants Jk with 
k > 4 and in a system with a monotonic dependence of 
Jfe on k is not very probable. 

An additional mechanism favoring the two-transition 
scenario may be related to quantum fluctuations, whose 
role in decreasing the energy of double walls may be more 
prominent than in decreasing the energy of a diluted net- 
work due to a smaller size of moving objects in the former 
case. 

The numerical simulations of the triangular-lattice 
Ising antiferromagnet have been performed in Refs. 
Ill El El In particular, Glosli and Plischke have 
studied the system with first and second neighbor in- 
teractions satisfying J2/J1 = 0.1, Rastelli et al. |2a 
- with J2/J1 =0.1, 0.5, 1, whereas Novikov et al. [25j 
have assumed that the interaction decays with the dis- 
tance exponentially. The results of these groups give ev- 
idence for the existence of a single first-order transition, 
which is consistent with our conclusions. In the simula- 
tions of Takagi and Mekata 27] the disordering of the 
striped phase has been investigated in the system with 
J2/J1 — 0.2 or 0.5 and J3/J1 = —0.75, but these authors 
make no conclusions about the type of the single phase 
transition which they observe. 

In addition to more traditional applications mentioned 
in the Introduction, the considered version of the Ising 
model can be used for the description of a triangular ar- 
ray of quantum dots at half-filling [2!| and (at sufficiently 
low temperatures) of a Josephson junction array with the 
dice lattice geometry and one-third of the flux quanta per 



plaquette [3^| . In the latter case the role of Jk with k > 1 
is played by the magnetic interactions of currents in the 
array [3I |3!| . The results of this work may also be of 
help for understanding the nature of phase transition(s) 
in the fully frustrated XY model on a honeycomb lattice 
03 l4l| , in which the fluctuation-induced vortex pattern 
[41| has the same structure as in Fig. ^ 

This work has been supported in part by the Program 
"Quantum Macrophysics" of the Russian Academy of 
Sciences and by the Program "Scientific Schools of the 
Russian Federation" (grant No. 1715.2003.2). 

APPENDIX A 

In the limit J\ — ► 00 it is convenient to split the Hamil- 
tonian into two terms, 

H = H + {H-H ), 

the first of which, Hq, includes only the infinite inter- 
action of nearest neighbors and restricts the summa- 
tion in the partition function to the states in which the 
constraint J2J) is satisfied on all triangular plaquettes, 
whereas the second term, H — Hq, includes all other in- 
teractions. It is well known that the application of a 
variational procedure allows one to use such a splitting 
to demonstrate that the free energy of the system, F, is 
bounded from above by 

F+ = F Q + (H-H )o, (Al) 

where Fq is the free energy of the system whose Hamilto- 
nian is equal to Hq, whereas the angular brackets denote 
the average calculated with the help of H . 
In our case in the thermodynamic limit 

F = -NiJx+Tso), 
(H-H )o = 3N( J252 + J353) , 

where N is the total number of sites, 

s re 0.323066 

is the residual entropy 0, whereas gk = (<7i<7j)o is the 
correlation function of the variables a on the sites i and 
j which are the kth neighbors of each other, calculated 
for the system with only nearest neighbor interaction at 
T = 0. According to Stephenson Q, 

12 13 

.92 - 7: + —/=- , 53 — 7: j ' 

9 V37r 9 TT 2 

The comparison of F + (T) with the free energy of a 
completely frozen ground state, which, naturally, coin- 
cides with its energy, 

E Q = -N{ Ji + J 2 - 3 J 3 ) , 

allows one to conclude that the temperature T c , at which 
a phase transition from a completely frozen state to some 
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phase with more developed fluctuations takes place, can- 
not be larger than 

T+{.h 1 J 3 )=C 2 J2-C 3 J 3 , (A2) 

where 

C 2 = (l + 3.g 2 )/s ~7.54, (A3) 

^3 = 3(1 -sa)/* « 11-08. (A4) 

Naturally, this approach does not allow to distinguish if 
the phase transition at T — T c < T+ is a direct tran- 
sition into the disordered phase or a transition to the 



intermediate phase in the framework of a two-transition 
scenario. 

For J 3 = the value of T+ following from Eqs. (|A2|) 
and (|A3(I is rather close to T 2 , the temperature of the 
spontaneous formation of double domain walls given by 
Eq. J7J, whereas in the limit of i?dw (that is 
^3 — * ^2/2) one gets 

T+ « 2.00 J 2 , 

which is compatible with an estimate for the temperature 
of the phase transition given by Eqs. JHJ and 1)11(1. 
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